Influence of supercoiling on the disruption of dsDNA 
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We propose that supercoiling energizes double-stranded DNA (dsDNA) so as to facilitate thermal 
fluctuations to an unzipped state. We support this with a model of two elastic rods coupled via 
forces that represent base pair interactions. Supercoiling is shown to lead to a spatially localized 
higher energy state in a small region of dsDNA consisting of a few base pairs. This causes the 
distance between specific base pairs to be extended, enhancing the thermal probability for their 
disruption. Our theory permits the development of an analogy between this unzipping transition 
and a second order phase transition, for which the possibility of a new set of critical exponents is 
identified. 



Unzipping of double-stranded DNA (dsDNA) as a pre- 
lude to transcription is a basic process of lifc^. There 
are presently two mechanisms that are considered for 
modeling this phenomenon. One is the Peyrard-Bishop- 
Dauboix model, in which localization of energy via a 
Fermi-Pasta-Ulam scenario of coupled non-linear oscil- 
lators causes thermal bubbles to form spontaneously at 
specific locations along a dsDNA chain^. The other 
mechanism is modeled by representing proteins, which 
cause unzipping, as an external forc o^i^i^ . This latter 
class of models envisions dsDNA as composed of cou- 
pled flexible Gaussian chains. These models describe 
the unzipping transition as a second order phase tran- 
sition. In realityS. the persistence length of dsDNA is 
approximately 45A. It may be a reasonable approxi- 
mation to treat a long semi-flexible chain as an effec- 
tive Gaussian chain with the monomer length substituted 
by the persistence length in order to investigate global 
characteristics^. But a realistic attempt to treat the un- 
zipping process through the application of an external 
force to a small region of the dsDNA must necessarily 
take into account the semi-flexible nature of the system. 

There is an additional outstanding issue which compli- 
cates our understanding of the unzipping of dsDNA. It 
arises through the observation that one can estimate^ the 
base-pair force as ^ O{100pN), while the experimentally 
observed^ minimum (external) disruptive force is only of 
the order of ~ O{10pN). Thus there would appear to be 
an intrinsic source of energy that allows the dsDNA to 
overcome this mismatch. What subtle mechanism comes 
into play during base-pair separation to explain this mis- 
match? 

Supercoiling is an ubiquitous feature of semiflexible 
rodsiS. One might conjecture that supercoiling plays 
a role in the packing of dsDNA. But its role in the 
disruption of base-pair bonds has been paid only scant 
attention^^. In this Letter, we make the qualitative case 
that supercoiling energy can have an appreciable effect 
on this phenomenon, by providing extra energy to free 
energy barrier to disruption. 

Benhamiiii^ was the first to suggest a phenomenolog- 
ical model to describe the influence of supercoiling on 
base-pair disruption. He obtained the parameters for his 



model by calibrating it to one set of dsDNA data, then 
verifled that the same choice of parameters worked for 
another experiment. Nelson et al."^ also referred to super- 
coiling in passing. Hennig and Archilla* have studied nu- 
merically an extended Peyrard-Bishop-Dauboix model to 
account for stiffness in dsDNA. They focused numerically 
on the effect of an external force applied to dsDNA char- 
acterized by varying degrees of overtwist. Their model 
is similar to the Poland-Scheraga model with torsion, 
which considers dsDNA as a semi-rigid one-dimensional 
ladder—. 

The model developed in this Letter provides a deeper 
and more detailed, analytical insight into the mechanism 
of base pair disruption in dsDNA. Marko and SiggiaiiS 
described the supercoiling of DNA by appealing to an 
evocative image of an over-twisted shoelace. The analogy 
that best describes the results of our model is that ds- 
DNA is like a ladder made out of bamboo. If the ladder 
is twisted beyond its normal, equilibrium planar state, 
one can imagine that it would be easier to break the 
rungs. We assert that this is precisely what happens in 
dsDNA. Supercoiling, or over-twisting, energizes the ds- 
DNA, making it more likely that thermal fluctuations 
will cause a given base-pair to be disrupted. Our theory 
provides qualitiative insight into the experimental obser- 
vations that superhelicity can make pBR322 DNA and 
E. coli duplex unwinding elements (DUE) susceptible to 
strand separation even if the over-twisting spans only a 
few base-pairsii. These are presumed to be examples of 
generic strand separation behavior. Wang et ali^ term 
this phenomenon stress induced duplex destabilization 
(SIDD). Benham.^-'^ and Wang et al^^ provide a valuable 
tool by assuming an Ising-type model that accounts for 
superhelicity, and with which they can obtain results con- 
sistent with experimental observations. 

We propose a model for dsDNA which is ab initio com- 
pared to BenhamiiiiS. It is based on generalizing the sin- 
gle elastica (conformations of a single, force-free elastic 
rod) model of Shi and Hearslii^ to coupled elastic strands. 
Shi and Hearst begin with a Kirchoff rod, and manipu- 
late the static equations of equilibrium into a non-linear 
Schrodinger equation! Here the amplitude of the wave- 
function represents the curvature of the rod, while the 
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argument of the phase factor is related to the torsion of 
the rod. We then consider on this basis two elastic rods, 
each exerting an attractive force on the other at pre-set 
points along their backbones to represent base pairs. Our 
approach thus has the flavor of a Salpeter-Bethe equa- 
tion for coupled quantum particles. We derive a non- 
linear model analogous to the Kronig-Penney model in 
condensed matter physics to understand the effect of su- 
percoiling on the coupling of the two coupled elastic rods. 
The advantage of extending Shi and Hearst's formalism is 
that it allows us to use well-established methods for solv- 
ing quantum mechanical problems to gain insight into the 
influence of supercoiling on dsDNA unzipping. 

A straightforward description of our model follows 
through the prescription of an energy functional £. When 
extremized with respect to its two component fields (each 
representing a strand), the functional yields a pair of 
non-linear Schrodinger equations (NLSEs), coupled via 
a (normal) force: 

£ = j dsf{^/j{s)) 

= fiMs)) + fiMs)) - A/(^i(s), ^2(S)) 

f =r(s)(^-c+i|V(s)p)v(s) 
A/ = 2A/'(24)(s)|^i(s)| -f 2AA(i,2)(s)|V2(s)| 

AA(2,1)(S) =-AA(i,2)(s) 

■^(2,1) (s) ^^i(Ji6{s-Si) (1) 

where the superscript zero denotes the functional associ- 
ated with each elastica, V' = + ^^y^ with lo^ and LOy 
being the components of the curvature vector in the x and 
y directions respectiveljsi^. The phase of V' is related to 
the torsion r of the strand^^. We have considered a com- 
mon label s for the location along either chain. This is 
sufficient to study the effect of super-coiling of a double- 
helix. We note that two independent labels would be re- 
quired to address issues such as transcription, for which 
one strand is translated relative to the other. The gener- 
alization goes through in a straightforward fashion, but 
is an unnecessary complication for the current problem. 

The subscripts 1,2 refer to each of the two semiflexi- 
ble strands in our dsDNA, s refers to the location along 
the spine of the dsDNA; A/(2,i), A/'(i.2) refer to the at- 
tractive forces exerted by the strand labeled 2 on the 
strand labeled 1 and vice versa. These forces are located 
at discrete positions along the backbone of the dsDNA, 
as signified by the Dirac delta function b(x). Oi repre- 
sents the force due to the i*'' base-pair. Using the po- 
tential developed by Rasmussen et ali^, we find that this 
force ranges from zero at the minimum of the potential 
well, and rises to a maximum of about 160pA^ for AT 
(Adenosine-Thymine) base pairs at the inflection point 
of the potential before tailing off. The maximum force is 
approximately AlTpN for CG (Cytosine-Guanin) pairs. 
The well-depths in these two cases are of the order of 
room temperature. 



As described by Shi and Hearslji^, the parameter c is 
a constant of integration and is related to the energies 
of bending and twisting, c — £b — / 4:, where £b is a 
dimensionless energy due to bending, and Q is a dimen- 
sionless energy associated with torsion. We shall treat it 
as an eigenvalue to be calculated to satisfy certain bound- 
ary conditions. The energy £ is chosen in units of M £, 
where M ~ 2000pA^ is the characteristic modulus of the 
dsDNA, and £ ~ 3.4A is the average distance between 
base pairs. The distances are measured with respect to 
i. 

Equation n is based on a refined NLSE for an elastic 
rod (See Eqn (3.10) in the paper by Shi and Hearst^^) in 
which certain transformations have been applied to cast 
the original relevant equations into a simpler form. More- 
over, we have retained only the component of the force 
between strands which is normal to each strand, while 
ignoring the tangential force. This provides an adequate 
model for coupled rods. 

By casting the theory of Shi and Hearslii^ in terms 
of an energy functional £, much physics can be gleaned 
in a manner similar to the Ginzburg-Landau theoryi^. 
Consider only the polynomial part of the energy density 
for each strand. It can be written down schematically as 
follows (note that the curvature k — ±|i/)|): 

f polynomial ~ ~ 2A/'|k| (2) 




FIG. 1: Phase transition depicted schematically from 
a single-well to a double-well structure. 

Figure n shows that for c ~ JV — 0, (red curve) there 
is a sole minimum at zero curvature (k = 0). This 
says, in essence, that when there is no binding force act- 
ing on the system and hence no bending or twisting, 
the minimum energy state is naturally a straight line. 
This comes about due to the quartic, non-linear term in 
the expression which embodies the semiflexible nature of 
each strand. This is a qualitative difference between the 
present case (red curve) and the ones described by the 
green and blue curves discussed in the next paragraph. 

The green and blue curves in Fig. prefer to configu- 
rations which possess different amounts of super-coiling. 
In either case, there are two minima, corresponding to 
two equal and opposite values of the curvature. The at- 
tractive force dominates the behavior for small curva- 
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ture, while the quartic term is important at high cur- 
vature. This results in the double-well structure. The 
blue curve refers to c = —0.125 and N — 0.025 (approx- 
imately 50piV). The green curve refers to c = —0.15 and 
N = 0.0125 (approximately 25pN). The greater nega- 
tive value of c for the green curve indicates a relatively 
greater degree of super-coiling. The state at the origin 
K = can be thought of as a transition state (when the 
strand is linear). The disruption of a base-pair corre- 
sponds to the achievement of this transition state. The 
well-depths, or equivalently, the barrier height between 
the wells dictates the probability that the system can 
reach the zero-curvature state through thermal fluctua- 
tions at a given temperature. This is governed by the 
standard Boltzmann factor. The green curve (with more 
negative c and therefore higher torsion) displays a lower 
barrier height as compared to the blue curve, and there- 
fore is more susceptible to attaining the disrupted state 
(k = 0) via thermal fluctuations. 

The qualitative change in the structure of fpoiynomiai 
in the absence and presence of an attractive force is anal- 
ogous to the critical point behavior in the theory of sec- 
ond order phase transitions, as the system passes from 
a single minimum to a double-well structure. However, 
there is a major difference in the two theories, since the 
present case concerns an order parameter, viz., the cur- 
vature, which differs from the usual sort of order param- 
eter (concentration, magnetization, etc.) which appears 
in Ginzburg-Landau theories. Given that our model ex- 
hibits a transition reminiscent of a second order phase 
transition, one can obtain a certain amount of insight 
into the universality class it belongs to by appealing to 
a mean-field approximation. In this case, by extrem- 
izing Equation [3 one obtains a cubic equation for the 
order parameter k. One can examine its behavior as 
A/" — *■ and c — > 0, as it is the simultaneous vanishing 
of these two parameters which determines the transition 
point. In this sense each of these two parameters plays 
a role analogous to that of 6 = T — Tc, the tempera- 
ture deviation from the critical temperature, in the stan- 
dard Ginzburg-Landau theory. If A/" — faster than c, 
then one obtains the conventional exponent of 1/2 since 
K ^ c^l"^ — > 0. On the other hand, if c — *■ faster than 
A/", then k ^ M^^'^ — > and suggests an exponent of 
1/3. Note that though there is but one independent vari- 
able in the model, viz., the position along the backbone, 
one can obtain three-dimensional structures through the 
specification of the curvature and torsion of dsDNA. 

Rudnick and Bruinsma have utilized the one- 
dimensional Poland-Scheraga model with torsion to sug- 
gest that there is a difference between a second order 
transition and denaturation in DNA-iS. Benham^^ pro- 
posed an Ising-like model to describe base-pair disrup- 
tion. It should be possible to determine from his theory 
whether the unzipping transition can be described as a 
second order phase transition. However, his energy func- 
tional is highly non-linear in the order parameter, and 
therefore difhcult to analyze in a simple manner. We 



speculate that it may be possible to map his model onto 
a non-linear sigma model (a continuum version of the 
Ising model) before analyzing it further. This is beyond 
the scope of our paper. 

We mapped the conceptual idea of base-pair disruption 
discussed above directly to a double-helix. By writing 

IV'('S)p| — \J'4'p{s) "ippis), p = 1, 2, one can extremize with 

respect to the complex conjugate fields ip*, and obtain a 
coupled set of NLSEs: 

Vp,p' = l,2; p'^p{3) 

Note that upon setting ^"2(5) 0, the model reduces 
to that of Shi and Hearst (see Eqn.(3.10) in that paper). 
The coupling between the rods occurs through the force 
exerted on one rod by the other. 

As discussed in detail by Shi and Hearst, the amplitude 
of the wave-function is the curvature of the strand, and 
the argument of the phase factor is related to the torsion 
of the strand. As is well-known in differential geome- 
try, curvature and torsion determine any given curve in 
space. Our solution of the NLSEs Eqn. Ois based on the 
idea that the dsDNA is composed of two counter-twined 
helices. According to the theorem of Bertrand^^, a helix 
is defined by the fact that the ratio of its curvature to its 
torsion is constant everywhere. 

With this in mind, we point out that a plane wave 
ansatz solves the NLSEs equation |21 with A/'(2,i) = 
-^(1,2) = (the force-free case of uncoupled elastica) as 
follows: 



V'p(s) = f^pi'^p) exp(iTps) 

KpiTp) ^^2{C+T^) (4) 

hCp IS the curvature, the total torsion in the strand is 
^totai ^ rp + a/2, a = (2/A - 1)Q, A is Poisson's ratio, 
and Q is a measure of the torsional energyi^.. We see in 
this case that the ratio of the curvature to the torsion 
is a constant everywhere. This occurs because Tp and, 
therefore, the amplitude of the wave-function is indepen- 
dent of the variable s. As such we observe that a plane 
wave solution to our force-free NLSE is consistent with a 
helix. 

We envision the solution of the full NLSEs with non- 
zero attractive forces between them by connecting suc- 
cessive helical solutions across each delta-function force. 
This is achieved by integrating the NLSE for each strand 
over an infinitesimal interval across any given delta- 
function force. This yields a condition connecting the 
first derivative of ip on either side of the delta- function 
force. We thus get a complex transcendental equation 
which yields two conditions (real and imaginary parts 
of the condition considered separately) to connect the 
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amplitude and phase of ip on either side of the delta- 
function. These conditions can be solved in a recursive 
manner to get values for the torsion and the parameter 
c in successive segments of the chain. If the phase dif- 
ference between the solutions on either side of the delta- 
function force is given by 5„ at the n*'* node, then the 
solution is given by the following set of recursive rela- 
tions: 

tanOn — 

'^n—l — 1 
l^n = \/2(c„ -I- t2) 

T"nl(Cn-l +T2_i)sec((5„)2 

= — i^^^TTs::? 

Tn = T„_i -I- 5n (5) 

Knowing the parameters for the helix at one end of the 
chain as a boundary condition, one can march recursively 
down each strand. In this sense we have solved an initial- 
value problem. Given that the dsDNA is two counter- 
twined helices {k — > — k), the parameters for the other 
helix are immediately obtained. 

The cartesian co-ordinates of a helix representing each 
segment between two consecutive delta-function forces 
are given parametrically {6) = xi"' +a„cos(0), 

yn{9) = Vn^ +a„ sin(6'), z„(0) = Vn^+knO, where is the 
aziniuthal angle. Xn\yn\yn'^ serve to fix the starting 
point of the n*^ segment. The parameters a„ and fc„ are 
related to the curvature and torsion of the helix. fc„ is 
the pitch of the helix, and a„ is the radius of the cylinder 
around which the curve is wound. These parameters are 
related to the curvature and torsion as follows: 
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Armed with this formal solution, we now display two 
cases to illustrate the notion that adding extra twist to 
dsDNA facilitates the path to the disruption of base- 
pairs. The parameters chosen for our calculations are 
representative of dsDNA, but do not refer to any specific 
DNA sequence. This was done to display the concep- 
tual basis of the model and provide a base-line for more 
realistic calculations which must follow. 

In the first case (see Fig. the constants (boundary 
conditions) are c = —0.15, a = 0.2 and Tp = 0.425 for the 
first segment. This yields approximately 10 base-pairs 
per turn. The base-pair separation is approximately lOA 
as in normal DNA. Hence it was designated the ground 
state of a model dsDNA and the forces between the base- 
pairs were set to zero. This would be consistent with an 
effective potential representing base-pairs^ which has a 
minimum at the equilibrium configuration, and therefore 
a null force. 



In the second case (Fig. 0) we took Tp = 0.46, keeping 
the other boundary constants unchanged. Thus this is 
an over- twisted case. This causes the double- helix to be 
squeezed down such that the distance between the base- 
pairs increases as we proceed up the ladder. The force 
was taken to be approximately 20pN. Our estimate for 
the force is in qualitative concert with the potential of 
Kim et alS. Unzipping, or base-pair disruption in our 
model, is signaled by a tendency of the base-pairs to 
show increasing separation. It is important to note, how- 
ever, that the increase in separation, indicated in units 
of £ = 3.4A in the figure along each base-pair, is not 
monotonic. as this happens, the absolute value of c„ < 
increases in successive segments, signaling a lowering of 
the free energy barrier depicted schematically in Fig^ 
There is a clear tendency of the separation to increase, 
followed by a decrease. This indicates that supercoiling 
is not distributed uniformly, but rather, can be localized 
over a small region of dsDNA. This is consistent with the 
model of Rasmussen et al. in which vibrations get local- 
ized to generate thermal bubbles^. Indeed it is possible 
that supercoiling and vibrations and supercoiling act in 
concert to achieve base pair disruption. 

There exist examples of genes that are only expressed 
once the relevant section of DNA is in a super-coiled 
stateii. It should be possible within the framework of 
our model to predict whether such a sequence is an ini- 
tiation site for transcription. At the moment we can ex- 
plain qualitatively via this illustrative example, the ex- 
perimentally known fact that in pBR332 DNA and E. 
coli^, supercoiling enhances unzipping. The examples 
just discussed show that given a certain torsion and cur- 
vature, progression over a finite number of base-pairs can 
create a localized increase in base-pair separation. This 
will lead to a lower attractive force. Eventually, the prob- 
ability of a thermal fluctuation to a state of zero curva- 
ture (disruption) can occur with appreciable frequency, 
as conceptualized in FigQ] 

We have restricted attention for the moment to the 
description of only a short sequence of like base-pairs, 
as that is sufficient to display the localization of super- 
coiling. We are solving an initial-value problem in which 
the conditions at one end of the double-chain are uti- 
lized to march recursively down the sequence. Therefore 
adding an extension to any given sequence will not affect 
the preceding section of dsDNA. 

In summary, we have extended the single-strand theory 
of Shi and Hearslii^ to two interacting elastic rods con- 
nected to each other via attractive forces. Their language 
is thus developed to the point where we can provide not 
only an over-arching framework for the unzipping tran- 
sition as analogous to a second-order phase transitions, 
but also puts us in a position to predict the behavior of 
specific sequences in future work. The analytical solu- 
tion provided in this paper can be generalized to more 
complex forces acting between the strands by using nu- 
merical techniques. This can yield predictions for specific 
sequences occurring in a given dsDNA which can be ma- 




nipulated experimentally to change the supercoiling. 



FIG. 2: Double-helical structure displayed over five 
sections for ci = —0.15, ri — 0.425, a — 0.2. The distance 
between base-pairs is marked on each rung, in units 
of 3AA. We have approximately ten base-pairs per 
turn for this case and it represents the ground state 
of dsDNA. There is a light source at the top left, 
casting a shadow of the double-helix. 
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FIG. 3: Double-helical structure displayed over five 
sections for ci = —0.15, ri — 0.46, a = 0.2. There is more 
torsion in this case as compared to the previous fig- 
ure. Base-pairs are separated increasingly compared 
to Figl2j but in a localized fashion, like the thermal 
bubbles calculated by Rasmussen et al.^. The use of 
delta-function forces in the model permits kinks in 
the structure. 
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